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Abstract. Wc consider the variational formulation of the electric field integral equation (EFIE) 
on bounded polyhedral open or closed surfaces. We employ a conforming Galerkin discretization based 
on divp-conforming Raviart-Thomas boundary elements (BEM) of locally variable polynomial degree 
on shape-regular surface meshes. We establish asymptotic quasi-optimality of Galerkin solutions on 
sufficiently fine meshes or for sufficiently high polynomial degree. 
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1. Introduction. Let T be a piecewise flat (open or closed) orientable surface 
equipped with a conforming triangulation A4 = {K}, consisting of triangles. Through- 
out, uniform bounds on the shape-regularity of the cells will be tacitly taken for 
granted, see [25, Ch. 3, § 3.1]. For a fixed wave number k > 0, let 14 and Vfc stand for 
the scalar or vectorial single layer boundary integral operator on T for the Helmholtz 
operator —A — k 2 , see [19, Sect. 4.1] or [21, Sect. 5]. The bilinear form underlying the 
variational formulation of the electric field integral equation ("Rumscy's principle") 
reads (see [3,35], [19, Sect. 4.2] or [21, Sect. 7.2] for closed surfaces, "boundaries", 
and [16, Sect. 3] for open surfaces, "screens") 

a(u,v) := (Vfcdiv r u,divrv) r - k 2 (V fc u,v) r , (1.1) 

where, as discussed in [20], (•, -) r hints at a duality pairing, extending the -^(Im- 
pairing for tangential vector fields or functions on T. The variational problem is posed 
on the Hilbert space 

X = H- 1/2 (div r ,r) , (1.2) 

in the case of a boundary T = <9f2, SI C K 3 a Lipschitz polyhedron, or on 

X = {ueH~ 1/2 (div r ,r): (u, grad r w) + (div r u, w) =0V«eC°°(r)} (1.3) 

in the case of a screen T. The latter space can be understood as a space of divp- 
conforming tangential surface vector fields with vanishing in-plane normal component 
on the screen edge dT. We refer to [14,17, 18,20,21] for a definition and more informa- 
tion about H|| (divp, T) and other trace spaces. An in-depth discussion for screens 
is given in [16, Sect. 2]. In this article we adopt the notations of [17]. Further, the two 
situations of open and closed surfaces will be treated in parallel. 
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The EFIE can be recast as a linear variational problem for a(-, ■) on X: given a 
source functional / G X' it reads 



In order to ensure uniqueness of the solution to (|1.4[) , we make the following assump- 
tion [21, Sect. 7.1]. 

Assumption 1.1. In the case of a closed surface T = <9f2 we assume that k is 
different from a Dirichlet eigenvalue of the operator curl curl on £1. 

We opt for a natural boundary element (BE) Galerkin discretization based on 
conforming trial and test spaces Xjv C X. These are obtained by using H(divp,r)- 
conforming Raviart-Thomas spaces of variable local polynomial degrees px £ No, 
K € M, on the surface triangulation yV(, see Sect. [2]for a precise definition. 

Refinement of the BE spaces can be achieved by raising the local polynomial 
degrees pk (p-refinement) or reducing the sizes Kk of the cells of M. (/i-refinement). 
Thus, the proposed discretization qualifies as "/ip-boundary clement method (BEM)" . 

Roughly speaking, judicious /ip-refinement can be expected to offer exponential 
convergence of Galerkin solutions even when the exact solution lacks global smooth- 
ness [44] . hp-BEM approaches have been suggested for various boundary integral equa- 
tions [5,6,32,34,36,37] and are a natural idea for the EFIE as well. While convergence 
theory for ^-refinement is well established [3,4,16,42], the extension to /ip-refinement 
proved to be difficult, see [8,9] for partial results. 

This article fills the gap and proves the following convergence result that translates 
into an a priori error estimate in the natural "energy norm" provided that information 
about some smoothness of the solution u of (|1.4|) is available, cf. [42, Sect. 8] and [16, 
Sect. 4] for the /i-version, [9] for the p-version, and [8] for the hp-version with quasi- 
uniform meshes. 

Theorem 1.2. There is a constant Co > such that for any f £ X' and for 
arbitrary mesh-degree combination satisfying max ^ ^° ^ e Galerkin BE dis- 

cretization of (jl.4p admits a unique solution u^r g Xjy and the Galerkin hp- BEM 
converges quasi- optimally, i.e., 



Both constants Co and C may depend only on the geometry of T and the shape- 
regularity of the surface triangulation M. . 

We remark that the policy of the proof has a lot in common with recent proofs 
of discrete compactness for the p-version of edge elements [11, 12,41]. The main tools 
are the same, namely, the sophisticated mathematical inventions of regularizing lifting 
operators [27] (sec Sect.|4]below) and projection based interpolation operators [29,30] 
(see Sect. [5]). They pave the way for verifying the assumptions of an abstract theory of 
Galerkin approximations for non-coercive variational problems, see [15] and Sects. [3l 
|6] below. 

Building on these mighty foundations the present article cannot be and does not 
aspire to be self-contained, but will give detailed references to relevant literature. We 
refer to [7] for an earlier version of this paper whose analysis is based on a Hodge- 
decomposition of X which, due to regularity issues on non-smooth surfaces, requires 
a sophisticated projection based interpolation operator which is not needed in this 
paper. 



u e X : a(u, v) = /(v) Vv e X . 



(1.4) 



||u-ujv||x < C inf ||u-vjv|| x ■ 



(1.5) 
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In the sequel, generic constants, designated by C, Co, C\, etc., may depend only 
on the geometry of V and the shape-regularity of M. . They must not depend on cell 
sizes, local polynomial degrees, and any function. 

2. Boundary element spaces. Raviart-Thomas surface elements provide an 
affinc equivalent family of divr-conforming finite elements under the Piola transfor- 
mation, see [13, Sect. III. 3] and [43]. We write HT p (K) for the local Raviart-Thomas 
space of order p on the triangle K G M, and TH~ p ,o(K) C H (div, K) for the sub- 
space of local Raviart-Thomas vector fields with vanishing normal components on 
dK. Vector fields in the latter spaces will be identified with their extensions by zero 
onto the whole surface V. 

Given a polynomial degree distribution \j>k ■ Pk G No, K G M}, we define edge 
degrees according to the "maximum rule" 

p E := max{px : K G M, E C K} , E G £ , (2.1) 

where £ is the set of edges of M. As elaborated in [39, Sect. 3.4], Raviart-Thomas 
spaces can be split into local "edge contributions" and "cell contributions" . In detail, 
write ipi and ip2 for the piecewise linear, continuous "tent/hat functions" associated 
with the endpoints of some edge E C E. We introduce the edge space 7VT PE (E) as 

nr VB {E) :=span{curl r (V>?V>2)> a,/?GN, a + P = p E + l}. (2.2) 
These spaces HT PE (E) obviously satisfy 



u N £lZT PE (E) => suppuAr C [J{K : E C K} and div r uAr = . (2.3) 



Then, we define the boundary element space (oblivious of boundary conditions!) ac- 
cording to 

X w = TZT (M) + KT PE (E) + ]T KT pK>Q (K) . (2.4) 

EES KEM 

Here, the space TZTq(M) is the lowest order Raviart-Thomas BE space. Thanks to 
the maximum rule (|2.ip . the localized spaces ~Kn(K) := Xjv|^-, K G AA, fulfil 

KT PK (K) cXjv(I) Vif G M and • n E \ E = V PE (E) V-E G £ , (2.5) 

with n^ standing for an edge normal, and V p for the space of (multivariate) poly- 
nomials of degree < p, p G No- Now, we are in a position to introduce the hp-BEM 
Galerkin trial and test spaces: 

• we pick Xjy := Xjv for closed surfaces V, 

• we choose Xjv := X^r n X for screens, that is, in order to obtain X^r edge 
spaces for edges contained in dT are simply discarded as well as basis functions 
of TZTq(A4) associated with edges on dT. 

Note that for E C K, K G M., we may encounter p E > Pk, and, consequently, 

X N (K) <£ UT PK (K) c X N (K) ! (2.6) 
However, thanks to (|2.4| and (|2.3[) . we can take for granted 

divrXAr (K) = div r TZT PK (K) . (2.7) 
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From [13, §111.3, Prop. 3.2] we know that div r TZT p (K) = V V (K). Thus, by (|2~4]) . 
divr : X^r > Qn, where Qn C £ 2 (r) is the space of yVf-piecewise polynomials with 
degree px on every K £ M. 

By [39, Theorem 3.7], [2, Sect. 5.5], the Raviart-Thomas BE space X^r allows 
for a discrete scalar potential space Sn C C°(T) comprising continuous piecewise 
polynomial functions on M. such that Sjq\ E C V PE +i(E) for all E £ £ and the 
localized spaces Sn(K) = Sn\k, K £ A4, satisfy 

curl r 5jv(A') = X N (K) n H(div r 0, K) VK £ M . (2.8) 

Below, we make repeated use of transformation to the reference triangle ("unit 
triangle") K := convex{( ), ( ), (°)}, which is mapped to a generic K £ M. by the 
affine mapping $> K : K h-> K, *(£) := A K x + t K , A K £ R 3 < 2 , t K £ R 3 . Writing &* K 
for the associated co- variant pullback of tangential vector fields, we define spaces of 
functions K i— > R 2 , 

Xjv(^) := *1.X W (A) , (2.9) 

which, due to non-uniform polynomial degrees, may be different for different cells K. 
The relevant K will be clear from the context. 

Remark 2.1. For the sake of brevity we do not include Raviart-Thomas BE spaces 
on (uniformly shape regular) quadrilaterals and BDM-type BE spaces on triangles into 
our analysis. With slight alterations the approach of this paper covers these settings. 
Besides, curved elements can be treated by the usual mapping techniques. 

3. Splitting technique. Owing to the infinite-dimensional kernel of divr the 
bilinear form a from (|1.1D fails to be X-coercive, which massively compounds the 
difficulties of convergence analysis for Galerkin schemes, as discussed, e.g., in [21, 
Sect. 3] and [24]. An abstract theory for tackling a priori Galerkin error estimates for 
non-coercive variational problems like (|1.4[) was developed in [16,22] and, in particular, 
in [15, Sect. 3]. The latter article tells us that Thcorem l 1 . 21 will follow, once we establish 

(A) the existence of a stable direct splitting X = V © W such that a\yxv an d 
a\Wxw are both X-coercive and a\vxw an d a\wxv are both compact, 

(B) the existence of a corresponding decomposition = Vat+Wjv, Wjv C W, 
that is uniformly stable with respect to cell sizes and polynomial degree p, 

(C) the gap property 



v - Vn\\ x ^ r, h K . . 

sup mf — — -<Omax< . (3.1) 

HvaHIjc " K ypK + i 



We remark that the approximation property 



inf llu — vat II Y — > as max \ — — — > 0, (3.2) 

vjvSXn" iV|IX K ^p K + l ^ ' 

dubbed CAS in [15,23], is automatically satisfied for families of hp- Raviart-Thomas 
spaces on families of uniformly shape- regular meshes. 

Taking the cue from the numerical analysis of electromagnetic wave equations [39, 
Sect. 5], one might resort to the "L 2 (r)-orthogonal" Hodge-decomposition of X [18] in 
order to obtain a suitable splitting. For /i-version analysis this idea was successfully 
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applied in [16,19,22,42]. Yet, on non-smooth surfaces, smoothness of functions in 
the V-component may be poor, which causes substantial technical difficulties, cf. [7]. 
These are avoided when following the guideline that the analysis of boundary integral 
operators is often greatly facilitated by taking a detour via a volume domain, cf. [26]. 
This strategy yields decompositions with enhanced smoothness of the V-componcnt. 

More concretely, as in [40] and [15, Sect. 4.3.1], V and W arc constructed via a 
regularizing projection R : X i— > X. To define them, we intermittently visit volume 
domains abutting T. There the construction employs Ti^-regular vector potentials, 
see [39, Sect. 2.4] and [1, Sect. 3]: 

Lemma 3.1. For any bounded Lipschitz domain O C R 3 there are continuous 
mappings L : curlH(curl, O) i— > (H 1 ^)) 3 and L : curl H (curl, fl) i— > (H^fl)) 3 
such that curl L 3? = <& for all 3? G curl H(curl, O) and curl L 4" = 3? for all 
<& G curl H (curl, 0). 

The construction is different for boundaries and screens and yields different projec- 
tion operators R c and R G , respectively, fortunately sharing the same pivotal properties. 
To begin with, fix u G X. 

(I) Case of a closed surface V = dil [40, Sect. 7]: R c u := ((L*) x to) L, where 

^ i TTim\ -Aw = in , 

v ' gradiu-n = divpu on 1 . 

The fact that J E divpudS 1 = for each connected component £ of T guaran- 
tees <fr G curl H(curl, f2), and Lemma I3TT1 can be applied. 
(II) Case of a bounded open orientable Lipschitz surface T with boundary dT and 
unit normal vector field 

Assumption 3.2. There exist two bounded Lipschitz domains S7i , C R 3 
satisfying 

• Tii n n 2 = r, 

• := Oi U L U Q2 is a bounded Lipschitz domain with trivial topology, 

• Tcffli and T C dQ 2 - 

In words, T is the cut chopping the sphere- like O into two parts fii, Q2, see 
Figure 13.11 

The fact that J T divpudS 1 = makes it possible to define for i = 1, 2 

-Aw, = in Qi , 

Wi G i? 1 (ili) : gradwi • n = on dfij \ F , 

gradu>i • rir = divru on F, 



and then 



graduii in fli 
gradu>2 in Q2 



G H o (div0,fi) , 



because the normal component of <& is continuous across T. Hence, we can 
apply Lemma l3~T1 and set R Q u := ((Lo 3?) x nr)| r - 

By using elliptic lifting theorems, the continuity of L and Lo, and trace theorems 
we conclude, * = c, o: 

3C = C7(r)>0: ||R.u|| H v 8(r) <C||div r u|| ff _ 1/a(r) Vu G X , (3.3) 

1 /2 

where HY (r) C X is the rotated tangential trace space 
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Fig. 3.1. Screen V with attached domains Qi and Q2- Note the nontrivial topology of F and 
how it can be dealt with in the construction of Q. 

• of (i? 1 (fl)) 3 on T := dfl for closed surfaces [17,20], 

• of (if,} (SI)) 3 on the screen L, sec [17, Sect. 3.2]. 
Moreover, by construction, on T 

div r R» u = div r u Vu 6 X => R 2 = R* . (3.4) 

Now we are in a position to define 

V:= R»(X)cHf (r), W := (Id- R*)(X) by JP X n H-^divrO, T) . (3.5) 

In light of (|3.3[) . the continuous embedding H^ 2 (L) <^-> H _1 / 2 (divr,r) ensures sta- 
bility of the splitting. 

Note that the embedding V w L 2 (r) is compact by (|3.3|) and Rcllich's theorem. 
Thus, thanks to the (iJ 1 / 2 (r))'-cocrcivity (rcsp., (H^ 2 (r))'-cocrcivity) of the single 
layer boundary integral operator Vk (resp., Vfc), see [19, Prop. 2], [21, Lemma 8], [21, 
Lemma 7], and [16, Proof of Thm. 3.4], we infer the X-coercivity of a\ VyV an d a,\wxw- 

1 /2 

Again, appealing to the compact embedding V (Hj_ (r))', the compactness of 
a \Vxw an d a\ Wx y is immediate [21, Lemma 9]. This yields (A). 

Remark 3.3. Assumption \ 3.2\ is easily verified for piecewise smooth Lipschitz 
screens through extension in normal direction followed by patching holes by means 
of thick cutting surfaces in order to mend topological defects. Yet, to keep the paper 
focused, we will not elaborate on this, but prefer to retain Assumption 1 3. °A 

Remark 3.4. We recall from [17] that X is the natural tangential trace space of 
J?(curl, SI) for a closed surface T, and of ifo(curl, 0) for a screen T. 
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4. Smoothed Poincare lifting. For a domain Del 2 that is star-shaped with 
respect to a e D, the Poincare lifting 



(P a u)(x) := / ru(a + t(x — a))(x — a) dr 
Jo 



(4.1) 



provides a right inverse of the 2D divergence-operator div u := (|^- + ^-) for contin- 
uous functions: div P a u = u for all u e C°(D), see [29, Sect. 3]. In [27] M. Costabel 
and A. Mcintosh demonstrated how to mend the somewhat insufficient continuity 
properties of P a by local averaging: 

Assumption 4.1. D is star-shaped with respect to a ball B C D. 

Then define the smoothed Poincare lifting [27, Sect. 3] as 

(P u)(x) := f ^(o)(P a u){x) da , (4.2) 



where ip s C°°(R 2 ), supp(-0) C B, f B ip(x)dx = 1. We get the following powerful 
mapping properties from [27, Cor. 3.3]. 

Theorem 4.2. Under Assumption ^. 1\ the smoothed Poincare lifting P according 
to (|4.2[) provides a continuous operator P : H S (D) \— * (H S+1 (D)) 2 for any s € K and 
satisfies div P <p = ip for all ip G L 2 (D). 

A crucial property of the smoothed Poincare mapping is the preservation of the lo- 
cal boundary element spaces: the smoothed Poincare mapping P on the (star-shaped) 
reference triangle K fulfills, cf. [41, Sect. 3], [38], [27, Sect. 4.2], 



P(dWRT p {K)) c nr p {K) by 



P(divXjv(/0) CX N (K) 



(4.3) 



5. Projection based interpolation. Following [40] and [15, Sect. 4.3.1] again, 
local projection operators will be used to build a suitable splitting of Xjy. How- 
ever, p-refinement entails a more subtle approach that resorts to so-called commuting 
projection based interpolation operators, see [29-31], [39, Sect. 3.6], and [28] for a 
comprehensive exposition. Commuting projectors link different finite element spaces 
on A4, the spaces Sn and X^r in the current setting. Employing the relatively simple 
construction of [29] will be sufficient for our purposes and the following results from 
that article and from [10, 30] will be used: 

1. There are projection operators (with domains £>(■)) 

U x : V{U X ) C flj 1/2 (div r ,r) ^ ± N , (5.1) 

U s :V(Ils)cH\T)^S N , (5.2) 
Hq : L 2 (T) « Q N , (5.3) 

satisfying the commuting diagram properties [29, Prop. 3] 

curlr o IT5 = IIx ° curlr on T>(Hs) , (5-4) 
divr o IIx = IIq o divr on T>(Ylx) ■ (5-5) 

For an open surface T the interpolation operator complies with boundary 
conditions: 



n x (xnp(n A -))=x„ 



(5.6) 
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2. As typical for the finite clement interpolation operators, Tlx and II5 are 
strictly local in the sense that both these projectors can be obtained by patch- 
ing together purely local cell based projectors Uk.x and Hk,s, K g A4. 
This is because for any edge E of M. with in-plane normal % the traces 
IIju • ue\ e and IIs^l^ depend only on (u -tie) \e and <p|-g-, respectively. 
Furthermore, pullback commutes with local interpolation: 

n K.x *a- = *k ° n A -,x on V(TL K ,x) ■ (5.7) 

3. The projectors IT^s enjoy the approximation property (this follows from [10, 
Thm. 4.1] with a scaling argument) 

\<p - TL Kt s<p\ mw < C ^^-J Mh^(K) W> G HV\K) . (5.8) 

These facts can be used to establish a special projection error estimate for Tlx x, 
cf. [41, Sect. 5], [39, Lemma 4.6], [1, Sect. 4]. 

Lemma 5.1. With C > depending only on the shape-regularity of the triangle 
K G M there holds 



|u - U K<x n\\ L2{K) < CJ l|u|| H v»(x) , 



for all u G H X > 2 {K) with div r u e div r X N (K) 



Proof. Write P for the smoothed Poincare lifting (see Sect. [4j) on K. Fix K e M 
and pick u € H 1 / 2 ^) with divu G divXj\r(iT). This vector field is split according to 

u = P div u + (u — P div u) = P div u + curl2D f , (5-9) 
div-frce 

where curl2D denotes a rotated gradient and the existence of the scalar potential 
<p G G H (K) : ffripdx — 0} is a consequence of div(u — Pdivu) = 0, which 
follows from Theorem l4.2l Theorem l4 . 21 also supplies the continuity of P : H~ 1 ^ 2 (K) 1— > 
H^ 2 {K), which paves the way for estimating 

M ff 3/2 ( £) < C\curh D ^\ H i/2 { K) < C (||u|| H1/2(jf) + ||Pdivu|| H1/2(jf) ) 

< H u llHi/2 ( if) +C||divu|| H _ 1/2(if) < C|H| H1/2(if) , (5.10) 

where the first step is justified by interpolation between iJ 1 (if) and H 2 (K). Then, 
by the projector property of Tlf> x , imbedding (|4.3p . and the discrete nature of divu, 
there holds 

u - n^. x u = (Id — 11^ x ) Pdivu+(7d — Tift x ) curl 2 n <p 

« ' 

=0 

by El , , , , n v 
= curl2/3(ia- 11^ s )(p , 
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where we owe the last identity to the commuting diagram property (|5.4[) on K. This 
makes it possible to apply 



U - U K.X U 



1530. 

< C (PK + 1)~ 1/2 \^ H 3/2 { k) < C( P K + l)- 1/2 |u| Hl/2(/?) 



Here, switching to the semi-norm in H 1//2 (_ftT) can be justified by a fractional Bramblc- 
Hilbert lemma [33, Prop. 6.1]. Eventually, l|5.7p and a scaling argument take the 
estimate to the cell K. □ 

We have implicitly proved that the (p-dependent !) local projectors Hk,x ■ {u G 
Yl l l 2 (K) : divpu e divrXjy(A')} i— > X.n(K) are continuous with norm independent 
ofp- 

Remark 5.2. In fact, the projector Tlx is closely linked to the splitting (|2.4|) . 
From [29] and [28, Sect. 4-] we extract the particular form 

ii x = n + J2 u E(id - n ) + ]T n K (id - u E )(id - n ) , 

Ee£ K 

where TLq, 11^;, IIr- are suitable projection operators into TZ.Tq(A4) , TUT PE (E) , and 
TVT PK fi(K), respectively. 



6. Discrete splitting. Since divrR*Xjv = divrX^, (|5.ip confirms that the 
following definitions are valid for * = c, o: 

Vjv :=n x R*(Xjv) , W N :=(Id-IL x °R*)X N . (6.1) 

By the commuting diagram property (|5.5|) and (|3.4[) . we find 



divrllx R* ujv = Ilgdivr R* u^r = TIq divrujv = divrujy Vujy £ Xjv ■ 

6Qjv 

=> R* TL X R* = R* on X w . (6.2) 

Hence, Iljt R* : Xjy i— > Xjy is a projection, which confirms that Xjv = V^v + W^r. 
Stability, in the sense of 

||n J tR.uj V ||j f < C\\u N \\ x , (6.3) 

with C > depending on L and the shape-regularity of M. only, is another consequence 
of Lemma 15.11 together with divrllx R* ujy = divrujy. This latter property also 
implies Wj, C W = X n H" 1 / 2 (div r 0, L). This verifies assumption (B) from Sect. [3] 

It remains to establish (C), the gap property (|3.1|) . which will be an immediate 
consequence of the following lemma. 

Lemma 6.1. There is a constant C > depending only on the geometry ofT and 
the shape- regularity of Ai, such that for * = c,o 



\(Id - Tlx) R* UjvIIx ^ C max J K \\u N \\ x Vu N e X N 

K V Pk + 1 
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Proof. By construction, we know that divp R* u^r = divpu^v, which permits us to 
apply the estimate of Lemma T5. II to R» Ujv|rt, K G M: 

\\(Id - Tl K ,x) R* u n\\ 2 L 2 (k - ) < C pK ^_ i l R * u Jv|hV2(k) • 



Patching together the local projectors and using sub-additivity of the |-| H i/2-semi- 

If 



1 /2 

norm, we arrive at (we remind that R* G Hj_ (r)) 



|(/d-n x )R* ujv|| L2(r) < CmaxW hjc [R*Ujv1 h 



l/2 (r) 



JO / /j 

< o max , 



K 



K V PK 



- ||div r uj V || J j-v2 (r) 



Since divp((J<i — IIx) R* ujv) = 0, this is sufficient for the assertion of the lemma. 

We point out that when we deal with an open surface T, we recall that R un G X 
is guaranteed by the construction of Sect. El For a continuous tangential vector field 
u G X that is smooth on the faces of T, the constraint in (|1.3p implies vanishing 
in-plane normal components on dT. By locality of Tlx, this will carry over to IljfU, 
which means IIxu G Xat, cf. (|5.6p . Further, we know from Sect. [5J that R Q is 
in the domain of 11^. A simple density argument then confirms Tlx RoUjv G Xjv, 
without adjusting the interpolation operator 11^, and the above proof carries over 
unaltered. □ 

The gap property (|3.ip now immediately follows from the estimate of Lcmma l6.lt 



. f l|v- v^vIIa- ||R* vjv - Vjvjjx 
sup mf — j| < sup j| 

vjvevjv vev !l v Jv|lx vjvevjy II v amIx 

EHJ || R* vtv — IIx R* vjv || y 
= sup 1| ji 



Lemma|6lT] / fij^ 

< C max n 



K y pk + 1 
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